function diff = q_0_fun( m, param)

    if param.K==0 
        diff = exp((1-param.alpha).*2./param.sigma^2.*int_delta_OJS( m, param));
    else 
        diff = exp((1-param.alpha).*2./param.sigma^2.*int_delta_fun( m, param));
    end

end